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ABSTRACT 



o 

^ • We present measurements of the mean dense core hfetimes in numerical sim- 

c^^. ulations of magnetically supercritical, turbulent, isothermal molecular clouds 
(MCs), in order to compare with observational determinations. The mean 

^ . "prestellar" lifetimes are given as a function of the mean density within the cores, 

^ . which in turn is determined by the density threshold nthr used to define them, 

l/^ . The mean lifetimes are consistent with observationally reported values, ranging 



from a few to several free-fall times. We also present estimates of the fraction of 
cores in the "prestellar", "stellar", and "failed" stages as a function of nthr- Failed 



. cores are defined as those that do not manage to collapse, but rather re-disperse 

back into the environment. Due to resolution limitations, the number ratios are 



measured indirectly in the simulations, as either lifetime ratios (for the prestel- 
lar cores), or as time- weighted mass ratios (for the failed cores). Our approach 
contains one free parameter, the lifetime of a protostellar object Tyso (Class + 
Class I stages), which is outside the realm of the simulations. Assuming a value 
Tyso = 0.46 Myr, we obtain number ratios of starless to stellar cores ranging from 
4-5 at nthr = 1-5 x 10^ cm~^ to ~ 1 at nthr = 1-2 x 10^ cm~^, again in good 
agreement with observational determinations. We also find that the failed cores 
are generally difficult to detect, although the mass in these cores is comparable to 
that in stellar cores at nthr = 1-5 x 10^ cm~^. At nthr = 1-2 x 10^ cm~^ the mass 
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in failed cores is negligible, in agreement with recent observational suggestions 
that at the latter densities the cores are in general gravitationally dominated. 
We conclude by noting that the timescale for core contraction and collapse is 
virtually the same in the subcritical, ambipolar diffusion-mediated model of star 
formation, in the model of star formation in turbulent supercritical clouds, and 
in a model intermediate between the previous two, suggesting a convergence of 
the models at least at the level of the core lifetimes, for currently accepted values 
of the clouds' magnetic criticality. 

Subject headings: ISM: clouds — ISM: evolution — stars: formation — stars: 
pre-main sequence — turbulence 



Introduction 



The evolution of the dense cores within molecular clouds (MCs) on their route to forming 
stars, and the duration of the various stages of contraction and collapse in particular, are key 
ingredients in our understanding of the star formation process. The process of contraction 
and final collapse of a core is generally divided in the so-called "prestellar" phase, which 
includes the span between the time when the core is first detected and when a protostellar 
object forms in its deepest regions, and the "stellar" phase, which includes the Class and 
Class I protostellar stages (from the appearance of the protostellar object to the clearance 
of the surrounding material) . 

Observationally, several works have estimated the durat ion of the prestellar stage through 



either statistical or chemical methods (see the discussion in Ballesteros-Paredes et al. 2007 



and references therein). The former rely on measuring the number ratio of prestellar to 
stella r cores, and assuming that this ratio equ a ls the ratio of durations of these stages 



Beichman et al. I 1986 



Lee fc Mvers I Il999l : lOnishi et al. I l2002l : iHatchell et al. 1 12007 



J0rgensen et al. 1 120071 ). A summary of prestellar durations, o r "lif etimes", from several 



studies, has been recently presented by IWard-Thompson et al. I (120071 ). 



Theoretically, there are two main competing models that describe the evolution of the 
cloud cores. The so-called "standard model" o f quas i -static contractio n mediated by am- 



bipolar diffusion (AD; e.g., iMouschovias I Il976l . Il99ll : IShu et al. 1 119871 ) was the dominant 



paradigm until recently. It postulated that the cores in low-mass star-forming regions evolved 
quasi-statically in magnetically subcritical clouds, contracting gravitationally at the rate al- 
lowed by AD, which causes a redistribution of the magnetic flux, until the central parts 
of the cores become magnetically supercritical and proceed to dynamical collapse. The 



3 



process was originally thought to be slow, because the clouds were thought to be strongly 
magnetically subcritical, in which case the AD timescale tad is roughly one order of magni- 
tude larger than the free- fall time. However, it was later recognized that near -magnetically 



critical initi al conditions (IFied. 



conditions ( Fatuzzo fc Adams 



er 



fc Afouschoy"iaslll993l : I ciolek fc BasulboOlh or turbulent 



2OO2I: iHeitsch et al. 



20041 ) cause the AD t imescale to be 



come significantly shorter, only a few times longer than the free-fall time ( ICiolek &: Basu 



2OOII : IVazquez-Semadeni et al. Il2005al . hereafter Paper I). This realization was motivated 
in part by o bservational studies suggesting core life times not much longer than their free-fall 



times (e.g., iLee fc Myers Ill999l : iJijina et al. Ill999l ). A recent numerical study of core life- 



times dominated by AD but with turbulent initial conditions has been presented recently by 
Nakamura fc Li I fl2005h . finding core lifetimes in the range 1.5-10 times the free-fall time. 



O n the other hand, the so called "turbulent" model of star formation fe.g..lBallesteros-Paredes et al 



1999al : IVazquez-Semadeni et al. II2OOOI : iMac Low fc Klessen 112004 : iBallesteros-Paredes et al. 



20071 ) takes into account the fact that MCs are turbulent, and attributes a dual character to 
the turbulence: on the one hand it provides support against generalized cloud collapse, but 
on the other it produces strong local density enhancements, which constitute the clumps and 
cores within the clouds. Thus, their formation timescales are of the order of the turbulent 
crossing ti me across the size scales of the fluid parcels that collide to produce the clump s 
and cores (IBallesteros-Paredes et al. I Il999bl : lElmegreen I I2OOOI : iPadoan fc Nordlund 1 12002| ) . 
However, Paper I showed that the lifetimes of randomly selected cores at densities compa- 
rable to those of observed ammonia cores are nevertheless of the order of a few times the 
free-fall time of the cores, because their evolution while they are detectable includes part of 
the assembly process and then the collapsing stage. Note that this is contrary to frequent 
perception s that the turbulent model im plies core lifetimes of approximately one free-fall 
time (e.g.. IWard-Thompson et al. 1120071 ). 



Within the context of the turbulent model, it is natural to expect the existence of 
an additional class of MC clumps that end up redispersing rather than collapsing; that is, 
turbulent density fluctuations that do not manage to become locally gravitationally unstable 
(Elmegreen 1993; Taylor et al. 1996; Paper I). We refer to these as "failed" clumps. Also, in 
this model, it is not crucial whether the clouds are magnetically subcritical or supercritical, 
since much of the "filtering" of mass that reaches the collapsed objects, and which regulates 
the star formation efficiency, is provided by the dual role of the turbulence. A stronger 
magnetic field simply appears to provide an extra reduction factor for the star formation 
efficiency, whose precise magnitud e depends on whether the turbulence in the clouds is driven 
or decaying (see the discussion in I Vazquez- Semadeni 1120071 . and references therein). 



In this paper we present a survey of the core lifetimes in a set of numerical simulations 
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of a turbulent, continuously driven, magnetically supercritical, isothermal MC, with the aim 
of obtaining statistically reliable data about them within this model, that can be compared 
to observational determinations. We also discuss the number ratio of "starless" to "stellar" 
(in the simulations, collapsed) cores. The plan of the paper is as follows: In §5] we describe 
the numerical simulations we analyzed; in ^we describe the method we used to measure the 
lifetimes and the number ratios of the various kinds of objects; in §4] we compare them with 
existing observational and theoretical work, and discuss their implications. In this section 
we also discuss the limitations and accomplishments of the present work. 



2. The Numerical Simulations 



We have chosen the moderately supercritical simulation (named M10J4/3.1) reported in 
Paper I for analysis of the core evolution. We refer to this run as "Rl" in the rest of the paper, 
and its selection was based on the fact that observations suggest that MCs are nearly critical 
or moderately supercr itical ( Crutcher 1999h. The comp utational domain contained 64 Jeans 



masses, or 1860 Mq (iVazquez-Semadeni et al. Il2005bl ). The mean number density is uq = 
500 cm~^, the simulation box size is 4 pc, and the (uniform) temperature is T = 11.4 K, so 
that the sound speed is Cg = 0.2 km s~^. The mean magnetic field was Bq = 14.5 jiG, giving 
a mass-to-magnetic-flux ratio = 2.8 times the critical value. These physical properties and 
conditions of the simulation make i t directly compara ble to the most massive "supercores" 
in the Perseus cloud, as defined by iKirk et al. I (120061 ). Indeed, supercore No. 2 (IC 348) is 
reported by those authors to have a mass ~ 1940 Mq, and, at a distance of 250 pc, its size 
is ~ 5 pc, very close to the corresponding parameters of our simulation. Supercore No. 7 
(NGC 1333) is reported to be somewhat smaller and less massive, at ~ 970 Mq. 

The turbulence in this simulation was continuously driven in order to maintain a tur- 
bulent Mach numb er ~ 10. For the turbulence random driver, we follow the method in 
Stone et al. I (119981 ) . i.e., we applied a random Fourier driver operating at the largest scales 



with the functional form P{k) oc A;® exp (— 8A;/A;pcak); where k = 27r// is the wavenumber 
corresponding to the scale /, /cpeak = 2(27r/L) is the wavenumber where the input power 
spectrum peaks, and L is the one-dimensional size of the computational box. As discussed 
in Paper I, this is motivated by the fact tha t turbulence is ubiquito us at all scales in molec- 



ular clouds (e.g., iLarsonl Il98ll : iBlitz I Il993l: iHever fc Brunt 



2004r). inclu ding starless ones 



such as the so-called Maddalena's cloud (IMaddalena &: Thaddeus Ill985l ). which suggests a 
universal origin and maintenance mechanisr n for the clouds' turbu l ence, and that turbulence 
is driven at the largest scales in the clouds (iHeyer fc Brunt II20041 . 120071 ). 



Two more simulations (called R2 and R3) were performed with the same physical pa- 
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rameters but different random seeds for the Fourier driver. R2 was performed at a resolution 
of 256 cells per dimension, and R3 was pe rformed at 512 ce lls per dimension. We used a 
total variation diminishing (TVD) scheme (IKim et al. Ill999l ). Runs Rl and R2 were inte- 
grated over 10 Myr, before self-gravity was turned on (5 turbulent crossing times). This is 
a standard procedure aimed at allowing the system to attain a well-developed, stationary 
turbulent state by the time gravity is turned on, thus preventing it from "capturing" features 
produced directly by the artificial random driver. Note that , however, one turbulent crossin g 
time is generally enough to achieve such a stationary state (IBallesteros-Paredes et al.ll2006l ). 
Thus, run R3 (the more expensive of these simulations in terms of computational time) was 
run only for one dynamical time before self-gravity was turned on. 

Once self-gravity was turned on, runs Rl, R2 and R3 were integrated over 13, 9.1, and 



population in the Solar Neighborhood 


(Bal 


esteros-Paredes et al. 


1999b: 


Hartmann et al. 


2001; 


Ballesteros-Paredes & Hartmann 


2007 


), molecular clouds should live not much longer 



than 5 Myr, since they are destroyed by stellar winds, bipolar outflows, SNe explotions, 
etc. For this reason, we only used the first 8 Myr of the evolution of runs Rl and R2 with 
self-gravity. 

Note that run Rl was used to perform the analysis over all types of cores (collapsed and 
failed). However, since only two collapse events were found in Rl, we used runs R2 and R3 
to verify that the timescales for formation and collapse obtained in Rl were consistent, and 
independent of the resolution. Run R2 produced two collapse events, while run R3 produced 
three. 



3. Measurement procedure 

3.1. General considerations 

Just as with observations, estimating the core lifetimes in the simulations is not a 
straightforward task. The numerical simulations have both advantages and disadvantages 
compared to the observations. First, in the simulations, the density, velocity and magnetic 
fields defined over the three spatial and one temporal dimensions are given. However, the 
cores, defined as we describe below, are elusive entities that in general do not preserve 
their identity over time; i.e., they do not involve the same mass at different times. After 
all, defining a core amounts to defining a certain region of space in what is really a fluid 
continuum. We choose to define a core as the set of connected grid cells surrounding a 
local density peak and having densities above a given density threshold rithr- This definition 
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mimics the emission observed when using tracers sensitive to different density reg imes. We 



have chosen not to use an algorithm such as CLUMPFIND flWilhams et al.lll994l ) because 
it can sometimes introduce artificial divisions across mostly continuous objects. Our chosen 
definition allows large-scale clumps to be defined at low nthr even if the region contains more 
than one local density maximum. 

In general, the cores are more poorly defined as lower values of nthr are considered, while 
instead this problem disappears as one considers higher nthr- Once defined, the cores must be 
followed over time. The animation of Rl associated to Fig. [T] in the electronic edition shows 
the evolution of the cores when nthr = 30 no = 1.5 x 10^ cm~'^, illustrating how the failed 
cores appear and disappear, and sometimes merge or split. It also illustrates a merger of 
collapsing cores (cf. §3.21) . The measured core properties (e.g., volume V, mass m, and mean 
density n) depend on nthr, and furthermore, at a fixed value of nthr, these properties vary in 
time for every core. Note in particular that the cores are not generally composed of the same 
fluid parcels throughout their history. Peak densities and their associated coordinates within 
a core do not depend on nthr, which only defines the time at which these two properties begin 
to be measured. 

The cores in the simulations can either collapse or rebound (no stable hydrostatic solu- 
tions are possible in the isothermal, magnetically supercritical regime; see Paper I). Thus, 
there exist two classes of starless cores: those that are on route to collapse, to which we refer 
as "prestellar" , and those that will eventually redisperse back to the surrounding medium, 
to which we refer as "failed" cores. In addition, we refer to cores that have already collapsed 
as "stellar" . With these considerations in mind, we see that the number ratio of starless to 
stellar cores is given by 

where iV^iess, -^pre, and A'^f are respectively the numbers of starless, prestellar, stellar and 
failed cores. 

Our goal in the present paper is to determine the lifetimes of the cores and their number 
ratios as a function of the threshold used to define them or, almost equivalently, as a function 
of their mean densities. The lifetimes that are reported in observational studies correspond to 
the duration of the prestellar stage, and are normally estimated by measuring the number of 
cores in the prestellar and stellar stages, and assuming that the ratio of their numbers equals 
the ratio of their durations. However, it was shown in Paper I that the ratio of the numbers 
is an upper limit to the ratio of durations if ther e exists a population of failed cores tha t 



is not accounted for. Furthermore, as discussed in iBallesteros-Paredes fc HartmannI (120071 ). 
potential source count incompleteness in several observational surveys also contributes to 
make this an upper limit. 
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In principle, one could simply measure the numbers of the various types of cores (failed, 
prestellar and stellar) in the simulations. However, a number of limitations prevent us from 
doing so. We have found that the collapsed, or stellar, cores in our simulations are clearly 
under-resolved. Their masses are typically between 50 and 100 M©, while the smallest and 
densest cores in, e.g., Perseus, have masses ~ IM©, although they are clustere d in groups of 



up to several tens of them (see, for example. Figs. 4 and 5 of lKirk et al. II2006I ). Specifically 



we have extracted the masses of the clustered su bmillimeter (submm) cores in supercore No. 



2 and supercore No. 7 from tables 1, 2 and 3 of iKirk et al. I (120061 ). We find a total mass in 
these clustered submm cores of ~ 7 Mq for supercore No. 2 and ~ 32 Mq for supercore No. 
7. In comparison, the masses at the time of peak density saturation of each of the collapsed 
objects in Rl at rithr = 1.2 x 10^ cm~^ (emulating the density of the submm cores) are 61 Mq 
and 62 Mq respectively. Thus, each of our cores with collapsed objects corresponds most 
directly to a cluster of protostellar cores. 

Another important issue is that the number of objects fluctuates statistically in time. 
This suggests that we should consider appropriately weighted time averages in order to obtain 
the most representative numbers. As mentioned, our analysis in Rl was performed over the 
first 8 Myr of the simulation after gravity was turned on, except for the failed cores at the 
lowest nthr; for which a sufficiently large statistical sample was obtained considering only 
the first 4 Myr (due to the large number of these objects). Analysis in R2 and R3 consisted 
only of collapse time measurements taken over 8 Myr and 5 Myr (the duration of R3) after 
gravity is turned on respectively. 

The time interval between outputs in the simulations was 0.04 Myr, and we used four 
different values of the density threshold: rithr = 1.5 x 10^, 3 x 10^, 6 x 10"^, and 1.2 x 10^ cm"^. 
Higher threshold densities were not used because cores with n > 256 uq = 1.28 x 10^ cm~ ^ 



are n ot well resolved in all the simulations according to the Jeans condition (ITruelove et al. 



19971 ). as explained in Paper I. Thresholds lower than 20 uq = 10^ cm""^ were not investigated 
either, because at such low densities the cores become impossible to follow in time due to 
their poorly defined identities (see, for example, the animations accompanying Paper I). In 
the next subsections we describe how we estimate the numbers and lifetimes of the various 
kinds of cores. 
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3.2. Collapsing cores 

3.2.1. Time scales 

When a core collapses in the simulations, in practice its peak density rip^ does not 
increase to infinity, but rather reaches a saturated value rigat ~ 2-3 x 10^ cm~^ at 256^ 
resolution or nsat ~ 1 xlO^ cm~^ in the 512^ simulation, beyond which the collapse cannot 
be followed further by the code. In this state, the mass of the collapsed object is spread out 
over a few pixels in each direction. This process is illustrated in Fig. [21 which shows the 
evolution of the peak densities of the objects that collapse in the three runs. Upper, middle 
and lower panels represent the peak density as a function of time for runs Rl, R2, and R3, 
respectively. The saturated density of a collapsed core in general oscillates in time. This is 
because the code contains no prescription for "capturing" the mass in the collapsed objects 
(the analogue of a "sink particle" in SPH codes, for example), and therefore this mass still 
interacts with its surroundings. Also, rigat increases slowly on average, because the collapsed 
object continues to slowly accrete mass. 

We define the prestellar lifetime of cL core, ^pre 5 the interval between the time when 
it is first detected at a given threshold n^hxi and the time when its peak density reaches 
the saturated value. This implicitly assumes that any further subfragmentation of the core 
occurs essentially simultaneously, so that all protostars form at roughly the same time, at 
least within the precision of our temporal resolution. This definition of the prestellar lifetime 
neglects the duration of the final stages of collapse (from the time when the density saturates 
in the simulation to the time at which actual protostellar densities would be reached). This 
approximation should not introduce a large error, since the free-fall time at the saturation 
densities is 

Tff = (37r/32G'p)^/2 _ 2 X 10^ yr, (2) 

which is significantly shorter than the temporal resolution of 4 x 10^ yr between snapshots 
in our simulation. In any case, this slight underestimate of the true prestellar lifetime acts 
in the opposite direction of any delay of the collapse introduced by numerical diffusion, and 
thus the two effects should tend to cancel each other out (see the discussion in §4.2p . 

The measured prestellar lifetime obviously depends on nthr, since this quantity defines 
the starting time of a core's observed "life". Empiricall y, we have found that the mea n 



density in the cores scales close to linearly with rithr (e.g., IVazquez-Semadeni et al. 1119971 ) 



and so this means that the measured prestellar lifetime should be a function of the mean 



density as well, in qualitative agreement with observational data (jWard-Thompson et al. 



2003). 

We ignore mergers of collapsing cores in the measurement of r^y-c- That is, if two cores 
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appear at times ti and t2 > ti respectively, and they merge at a later time into a new core 
that eventually collapses, we take the prestellar lifetime Tprc as the time elapsed between ti 
and the time of density saturation. We do this motivated by the fact that the peak densities 
ripk of cores that merge into a new core which collapses later are seen to rise rapidly even 
before the merging event, suggesting that these cores were already on their own independent 
route to collapse, rather than the latter being triggered by the merging event. For instance, 
the first two cores in Rl (upper panel of Fig. [21 solid and dashed lines respectively) are seen 
to appear at roughly the same time (t ~ 1.4 Myr), and to merge shortly thereafter. Thus, 
we consider them together as the single first collapse event in Rl. The third core is seen to 
have very similar physical properties to those of the merger, and we refer to its collapse as 
the second collapse event in Rl. In Fig. [3] we report the prestellar lifetimes measured for 
the collapse events in Rl (diamonds) , R2 (triangles) and R3 (squares) as a function of their 
mean densities (which in turn are ~ nthr)- 

It is important to remark that the long durations of the collapsed objects after they 
have collapsed seen in Fig. [2] does not imply that the stellar cores last that long, because 
no modeling of stellar energy feedback (e.g., winds, bipolar outflows, ionizing radiation) is 
included in the numerical model, and thus the cores cannot be dispersed after they produce 
collapsed objects. 



3.2.2. Number ratios 

As mentioned in §3.11 we cannot simply count the number of cores with collapsed 
objects because it is clear that at the resolutions we used their subfragmentation was not 
appropriately captured. Thus, instead, we proceed inversely to the standard observational 
procedure, and estimate the number ratio of prestellar to stellar cores simply as the ratio of 
the prestellar lifetime to the lifetime of the embedded protostellar phase: 

TV r — 

-''pre 'pre /qx 

where 7^ is the average of the measured prestellar lifetimes for the seven collapse events in 
the simulations, and Tyso is the duration of the protostellar stage. We emphasize that eq. 
([3]) allows us to estimate the number ratio of prestellar to stellar cores, but not the actual 
numbers, and that this ratio is a function of the threshold density, since Tpre is. 

For US, '^yso IS cL free parameter because it is outside the realm of our simulations. 
As mentioned in §3.2.11 no modeling of the feedback from the protostars such as bipolar 
outflows, ionization radiation or winds is included, so there are no energy sources available 
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for dispersing the cores once they are formed. Thus, we have to consider it as a free parameter 
in our study, whose value is to be set on the basis of external information. 

Observationally, Tyso is normally estimated by o btaining the ratio of the number of 



embedded objects to that of the T Tauri stars (e.g., IWilking et al. 



1990l : iGreene et al. Ill994l : iKenyon &: Hartmann Ill995l : iHatche 
reported in the literature are 



1 et al. 



19891: iKenvon et al. 



20071). Values fo r r, 



Ward-Thompson et al. 112007 



0.6 ± 0.3 Myr (see review by 
and references therein), the spread being due to both observational constraints and real 
differences among the observed regions. 



We take Tyso = 0.46 Myr, the mean value obtained for Perseus by lHatchell et al. I ( 120071 ). 



We base this selection on two considerations: First, their estimate of Tyso includes corrections 
for incompleteness of surveys toward very low masses and multiplicity. Second, we have 
found that o ur simulations are reasonably resemblant of a massive "supercore" in Perseus, 
as defined bv lKirk et all J2006h (see g2]and ^3J\f . 



3.3. Failed cores 



Failed cores in our simulations are less dense than collapsing cores, and therefore are not 
as strongly affected by the limited resolution of our simulations. Nevertheless, we consider 
it safer to estimate their number ratio with respect to the stellar cores through the ratio of 
masses in these two kinds of objects. This estimator avoids the uncertainty associated with 
the possible failure to capture further subfragmentation of the cores. However, even this is 
still not completely straightforward because we need to take into account temporal effects. 

The failed cores "appear" at some time and "disappear" at a later time at a given 
density threshold nthr- Furthermore, their masses are continuously varying over this time 
interval, first increasing, reaching a maximum value, and then decreasing again. In order to 
obtain an estimator that is most representative of the number ratio expected to be found 
upon a single instantaneous observation, we take the ratio of the temporal averages of the 
masses of the failed cores to those of the stellar cores, weighted by the ratio of their typical 
durations 



ry,oE,K)' 

where the summations over i and j refer to the failed and stellar cores respectively, (mj) 
and Ti are respectively the time-averaged mass and the duration of the i-th failed core, and 
Tyso is the duration of the protostellar phase discussed in §3.2[ For the stellar cores, we take 
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(mj) as the average mass over a time interval Tyso after the core has collapsed. 

Equation (jl]) assumes that actual failed cores have, on average, the same mass as stellar 
cores. However, it is reasonable to expect failed cores to be less massive in general than 
individual stellar cores, and in this case the right-hand side of eq. (jlj) is a lower limit for 
the failed-to-stellar core ratic0. Furthermore, eq. H] also assumes that no fragment of a 
failed core ever ends up collapsing, and no fragment of a stellar core ends up redispersing. 
This is a reasonable assumption, since subfragmentation of the cores is mostly a gravitational 
phenomenon, which occurs after the onset of collapse. Thus, collapsing cores are expected to 
fragment into further collapsing units, while failed cores are not expected to fragment because 
they do not engage into collapse. In fact, the failed cores in our simulations generally have 
very small, sub-solar masses (see §4.ip . 



3.4. The starless to stellar ratio 

Equations ([3]) and (jlj) give us the prescriptions to estimate the terms on the right-hand 
side of eq. ([1]). We can then estimate number ratios with respect to the total number of 
cores A^tot = A^pre + A^'f + A^* in terms of the ratios derived in eqs. ([I]), ([3]), and (jlj) as follows: 



N 

^ "pre 


r Artless /A^. 


A^tot 




A^f 








A^+less 
iVtot 


1 f 

^ V A^. 


A^. 
A^tot 


-| A^ *less 



pre 



A^, 
A^f 



(5) 
(6) 
(7) 



These fractions are shown in Fig. H] {left panel) as a function of nthr- The right panel of this 
figure shows the same ratios but normalized to the number of stellar cores. 



^See, however, the discussion in §4. II concerning the possibihty that failed cores are hard to detect even 
at low densities, an effect that would act in the opposite direction, lowering the observed failed to stellar 
core ratio. 
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Discussion 



4.1. Comparison with previous work 

In addition to the core hfetimes as a function of their mean densities, Fig. [3] shows 
two lines indicating the locus of t = rg and t = lOrff, with Tff d efined in eq. ([2]). This 



figure can be compared with Figure 2 of IWard-Thompson et al. I (120071 ). It is seen that 



within the uncertainties, the core lifetimes in our simulations are in good agreement with 
the observational values, being in general a few to several times the free-fall time. However, 
such a comparison should be taken with care because in the determination of observational 
lifetimes, the assumed values for Tyso span one order of magnitude (from 10^ to 10^ yr). In 
any case, one can unambiguously conclude that the prestellar core lifetimes in the turbulent 



scenario are a few times Tff, and not only Itq as is often stated in literature (e.g.. iKirk et al. 
20051 : IWard-Thompson et al. IboO?! ). 



These results can be understood as a consequence of several factors. First, the cores 
take some time to be ass embled by the extern al turbulent compressions before they become 
gravitationally unstable (IGomez et al. 1120071 ). and even after that, some readjustments of 
the mass and some oscillati ons are expec ted to occur before gravitational collapse sets in. 
Moreover, as pointed out by iLarsonl (119691 ). the duration of the actual collapse is longer than 
the free-fall time because the thermal pressure gradient is never negligible. So, even though 
the cores in the simulation are formed dynamically, their lifetimes are fully consistent with 
observational estimates. 

The number ratio of prestellar to stellar cores in our simulations and the ratio of dura- 
tions of thes e stages at low value s of rithr are ~ 3 (Fig. HI right panel). The same value is 
obtained by iLee fc Myers I (119991 ). who studied an optically-selected sample of cores, char- 



acterized by relatively low mean densities, < 10^ cm~^. Note, however, that we obtain 
a number ratio of total starless (prest ellar + failed) to stellar cores of ~ 5, larger than 
the observation of iLee fc Myers I (119991 ). This suggests that their sample consisted mostly 
of prestellar cores, possibly because the failed cores tend to be smaller and therefore less 
easily detected. Indeed, Fig. [5] {left panel) shows a size histogram of the failed cores at 
'^thr = 30 rio = 1-5 X 10^ cm~^ for Rl, and it can be seen that vi rtually all of them ha ve sizes 
£ < 0.2 pc, with £ ~ 0.07 pc being a characteristic size; instead, iLee Sz Myers I (119991 ) report 
a mean size £ ~ 0.2 pc. Thus, real failed cores may be harder to detect observationally, and 
therefore not very common in observational surveys, even of low-density cores. For reference, 
a histogram of the failed core masses is shown in the right panel of Fig. [5] as well. 

At high thresholds, rithr = 1-2 x 10^ cm~^, the ratio Npj.c,/N^ ~ 1.4. This is consis- 
tent with the observed value of ~ 1 for this ratio in recent surveys of dense, submm cores 
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( iHatchell et al. 1120071 : 1 J0rgensen et al. 1120071 ). sensitive to n ~ 10^ cm ^. Thus, with a single 



assumed value of Tyso, our estimated number ratios are in good agreement with observations 
at both low and high densities. Moreover, the number of failed cores at the higher densities 
goes to zero, suggesting that indeed submm surveys sample essentially all gravitationally 
bound cores, thus tracing the final stages of core evolution prior to protostar formation. 

Our estimated core lifetimes are also in good agreement with the predictio ns from the 
AD-m ediated model a pplied to moderately sub critical clouds, either in analytical (ICiolek &: Basu 
200ll ) or in numerical (INakamura &: Li II2005I ) form. The latter authors report lifetimes 1.5- 
10 times the local free-fall time in the collapsing cores, while the lifetimes of our cores are 
~ 6 Tff (see Fig. [3]). The similarity between the mean lifetimes in both cases occurs because, 
on the one hand, tad is only a few times longer than Tg at the densities and observed values 
of the mass-to-magnetic flux ratio of molecular cloud cores while, on the other, the cores 
in our simulations are not necessarily in a direct route to collapse when first detected, but 
instead undergo a period of build up and readjustment, as mentioned above. Therefore, 
we conclude that both models give comparable predictions for the lifetimes of the cores in 
their observable stages, with AD not significantly delaying the final stages of contraction 
and collapse. 



4.2. Caveats, limitations, and error estimates 

The study presented in this paper has faced a number of difficulties comparable to those 
encountered by observational studies, and in this sense, its results are only suggestive, rather 
than conclusive. The main limitations were: 

1. The limited resolution of the simulations. This limitation prevented us from correctly 
following the fragmentation of collapsing cores, and forced us to estimate the number ratio 
of prestellar to stellar cores through the ratio of their durations (rprc/ryso) rather than through 
direct counting, with Tyso being a free parameter of our study. Similarly, the number ratio 
of failed to stellar cores had to be estimated in terms of the ratio of their masses, again due 
to the inability to resolve individual stellar cores. 

The limited resolution also restricted the number of collapse events to only ~ 2 in each 
of our runs, although each collapse event clearly corresponds to the formation of a stellar 
cluster (i.e., to many collapse events leading to individual stars), since the masses of the 
collapsed objects are 50-100 Mq. Unfortunately, we estimate that, in order to adequately 
resolve these objects into solar-mass- like cores, an increase of at least a factor of 50^/"^- 
100^/^, or ~ 4 in resolution would be required (i.e., resolutions > 1024^), or the usage 
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of a Lagrangian-type of numerical algorithm (SPH or adaptive-mesh- refinement), in both 
magnetic and non-magnetic regimes. We hope to perform such a study in the near future. 

The limited resolution might also possibly affect the collapse times measured in the 
simulation (see §3.2.11) . both because of the saturation values of the density reached by the 
collapsed objects as well as because of the possible slowing down of the collapse due to 
numerical viscosity. To test for this, Figures [6] and [7] illustrate the dependence of collapse on 
resolution. Figure E] shows the evolution of the maximum density rimax in four simulations 
with the same global parameters (turbulent Mach number ~ 10, Jeans number J = 4) at 
resolutions 64^ [red line), 128^ {green line), 256^ {dark-blue line), and 512^ {clear-blue line). 
The collapse in each simulation consists of the sharp transition from a low typical level of 
the maximum density {rimax ^ 100 uq = 5 x 10^ cm~^) to a high one {rimax ^ 10^ tiq = 
5 X 10^ cm~^). We define the "collapse time" of the simulations as the duration of the sharp 
rise of nmax over more than one order of magnitude separating the two levels. This is well 
defined for resolutions equal or larger than 128^. For economy reasons, the simulation at 
512^ (i.e., R3) is only run for a short time compared to the other simulations. In all cases, 
there is seen to be a rebound from the maximum. 

Figure [7] summarizes the collapse times as a function of resolution. For the 128'^ sim- 
ulation the collapse time corresponds directly to the measured duration of the rise of nmax 
for the first collapse, and the error bar denotes the full width of the rebound around the 
saturated value. For the 256^ (Rl) and 512^ (R3) simulations, the collapse times correspond 
to the average of the durations measured for the available collapse events. The error bars 
denote the range of values found. This figure shows that the collapse time appears to be 
independent of resolution for resolutions 128^ and above, meaning that our longer-than-free- 
fall collapse times are not an artifact of the resolution. If anything. Figure [7] shows that 
there is a slight trend of the collapse time to increase with resolution, probably suggesting 
that the effect of traversing a larger dynamic range in density at higher resolution dominates 
over any slowing down of the collapse by numerical viscosity at lower resolution. 

2. The simulations do not include ambipolar diffusion (AD). Even though the simulations 
were supercritical, in the presence of AD there would exist the possibility that some failed 
cores could be "captured" by AD if the main agent supporting them is the magnetic energy, 
since in this case AD could cause a redistribution of the magnetic field, leaving the centermost 
parts of the core with less support than in the case without AD. This, however, does not 
appear to be the case. For example, the longest-lasting failed core in Rl, with a lifetime of 
~ 1.9 Myr defined at rithr = 3 x 10^ cm~'^, is clearly supercritical (/i ~ 4.0±2.1 times critical, 
where the uncertainty comes from the estimation of its "radius" , since the core is far from 
spherical; see Paper I), indicating that its failure to collapse was due to thermal-|-turbulent 
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rather than magnetic support. 

In any case, a worst-case estimate of the error committed by the lack of AD can be made 
by assuming that all failed cores with lifetimes longer than a representative AD timescale tad 
proceed to collapse rather than rebound. A reasonable estimate for tad is ~ 1.5 Myr (see 
the Appendix in Paper I). Being even more restrictive, we recalculated the number ratios 
assuming that the failed cores in Rl with lifetimes r > 0.9 Myr are "captured" by AD (see 
Fig. [H]). In this case, we find that the estimated number ratios of failed cores decrease by a 
factor of ~ 2, but this does not significantly affect the starless to stellar core ratio, because 
the failed cores have very low masses. 



4.3. Conclusions 



In spite of its limitations, the results of the present study are nevertheless encouraging, 
as the prestellar lifetimes measured directly from the density evolution of the cores are in 
good agreement with observational determinations. The same is true for the number ratios 
of prestellar to stellar cores. This estimation depends on the free parameter Tysn , but a single 
value for it, chosen as the mean of the range reported by lHatchell et al. I (120071 ) for Perseus, 
produces number ratios that are in good agreement with observational determinations at 
both high and low values of the density threshold n^Ya- The present study thus suggests that 
the turbulent model of star formation is not inconsistent with observational determinations 
of core lifetimes and the number ratios of prestellar to stellar cores. 
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Fig. 1. — Cores defined at nthr = 30 uq = 1.5 x 10^ cm~^, and at t = 1.96 Myr for Rl. Tliis 
time marks tlie end of collapse for the first collapsing core (the biggest core seen near the 
center of the simulation box) in this run. At this threshold, failed cores become an important 
feature and many of them coexist at every temporal step. The two "spots" at opposite sides 
of the y-axis boundaries, and at {x, z) = (203, 197) correspond to just one failed core crossing 
the periodic boundaries. An animation of the simulation for 0.32 < t/Myr < 8 is presented 
in the electronic version of the Journal, with the frames separated by At = 0.04 Myr. The 
first 9 frames are omitted, since there are no structures above rithr at those times. The figure 
shown in the printed version of the Journal corresponds to frame #40 in the animation. 
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Fig. 2. — Peak density as function of time for collapsing cores in Rl {upper panel), R2 {middle 
panel) and R3 {lower panet) after they have exceeded ripk = 6 x 10^ cm~^. A merging event 
occurs among the first two collapsing cores of Rl {solid and dashed lines respectively), but 
their peak densities are seen to rise steeply before their merger (the time when the lines 
converge). A merger also occurs in R2, but well after the individual collapse of the cores. 
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Fig. 3. — Duration of the prestellar stage for the collapsing cores in the simulations as 
a function of their mean density at the time collapse starts being measured. Diamonds, 
triangles and squares respectively correspond to the Rl, R2 and R3 runs. The loci of Tpre = Tg 
and Tpre = 10 Tff is marked by the lower and upper long-dashed lines respectively. Also shown 
is the fit for the mean prestellar lifetime as a function of mean density {dash- dotted line), 
which is ~ 6 Tff. 
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Fig. 4. — Left panel: Number ratios of the stellar {dashed line, diamonds) , total starless {solid 
line, triangles), prestellar {dotted line, plus signs) and failed {dash-dotted line, crosses) cores. 
Right panel: Number ratios of the total starless, prestellar and failed cores (same lines and 
symbols as in the left panel), normalized to the number of stellar cores. 
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Fig. 5. — Left: Size histogram of the failed cores at nthr = 30 no = 1.5 x 10"^ cm ^ for Rl. 
Right: Mass histogram for the same cores. 



- 23 - 




Fig. 6. — Evolution of the maximum density rimax in four simulations with the same global 
parameters (turbulent Mach number ~ 10, Jeans number J = 4) at resolutions 64^ {red 
line), 128^ {green line), 256^ {dark-blue line), and 512'^ {light-blue line). The time axis is 
given in units of the sound crossing time {lower axis), equal to 20 Myr, and of the quantity 
Lj/cs {upper axis), which is close to the free-fall time. The density axis is in units of the 
mean density in simulations uq = 500 cm~^. The first collapse event in each run is denoted 
by the sharp transition from a low mean maximum density level to a high one. The "collapse 
time" of the simulation is defined as the duration of the sharp rise of nmax over more than 
one order of magnitude separating the two levels. This is well defined in all cases except at 
64^. 
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Fig. 7. — Collapse time versus resolution for the three highest-resolution runs shown in Fig. 
[6l For the 128^ simulation, this time corresponds to the measured duration of the rise of 
'^max for the first collapse, and the error bar denotes the full width of the rebound around 
the saturated value. For the 256^ and 512^ simulations, collapse times correspond to the 
average of the durations measured for the collapse events shown in Fig. [21 The error bars 
denote the range of values found. 
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Fig. 8. — Same as Fig. Hlbut assuming that all failed cores with lifetimes r > 0.9 Myr are 
"captured" by AD and led to gravitational collapse. The fraction of failed cores is seen to 
decrease, but their masses are so small that the fraction of stellar cores remain essentially 
unchanged. 



